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We derive the gauge-free Hamiltonian structure of an extended kinetic theory, for which the intrinsic spin of 
the particles is taken into account. Such a semi-classical theory can be of interest for describing, e.g., strongly 
. - . magnetized plasma systems. We find that it is possible to construct generalized brackets for the extended phase 

space, and discuss the implications of our findings. 
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' High energy density plasma physics has become a popular subject (see, e.g., fst] and references therein). In such systems, 
quantum mechanical effects, such as wave function dispersion and/or statistical effects, can become important (for a recent 
d experimental example, see 1 1 1]), and much of these plasmas can be rightly termed quantum plasmas. Much of the early literature 
on quantum plasmas has focused on condensed matter systems with a background lattice structure and the linear effects that 
follows (see , e.g.. (H Oll). However, recent developments shows a different direction, where the nonlinear aspects of such systems 
are in focus Il20ll46ll . Examples of recent results include quantum ion-accoustic waves |2], Jeans instabilities in quantum plasmas 
. ll43ll . trapping effects jiill , magnetization by photons i45ll and relativistic effects |25, 26]. Typically, the quantum hydrodynamic 
(~| equations are derived by starting from the Schrodinger equation and making the Madelung ansatz 1 19]. However, a method that 
more closely resembles the classical case is to use kinetic equations as a starting point (see Ref. 1 19] for a comparison between 
"i the different approaches). The field of quantum kinetic theory |0] in many ways started with the ambitions of Wigner, as 
^ presented in Ref. |50], to bridge the gap between classical Liouville theory and statistical quantum dynamics 1 12, 36, 49]. Thus, 
^ the development of quantum kinetic theory was partly due to an interest in obtaining a better understanding of the quantum- 
,— — I to-classical transition lf52ll . However, another important aspect of quantum kinetic theory is as a computational tool for, e.g., 
quantum plasmas fi^, condensed matter systems lfT4ll22ll . and, in general, quantum systems out of equilibrium |'4l'], and in that 
respect shares many commonalities with quantum optics flst]. As shown in [5, 6, 8, 17, 21], spin is such an effect, the one of 
. ^ particular interest here. 

When developing new models it can be important to show that they are Hamiltonian, since all of the most important models 
of physics have this property when phenomenological or other dissipation is neglected. Matter models in terms of Eulerian 
or spatial variables possess noncanonical Hamiltonian form, i.e. they are Hamiltonian but the conventional variables are not 
a canonically conjugate set and consequently the Poisson bracket possesses noncanonical form yet retains its Lie algebraic 
properties of antisymmetry, bilinearity, and the Jacobi identity. (See e.g. |28-30, 33]. Also see lfll l48ll for recent work applicable 
to plasmas.) 

For the kinetic theory of interest here we show it has a Hamiltonian structure that is a generalization of that given in ll24ll27ll28ll . 
(See also lHI].) We present the noncanonical Poisson bracket, prove directly that it satisfies the Jacobi identity, find Casimir 
^yl^ ' invariants for the theory, and present an energy-like theorem that demonstrates that all equilibria with monotonically decreasing 
' distributions are stable. 

We consider the nonrelativistic spin Maxwell- Vlasov equation for /(x, v,s,f), an electron phase space density: 
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where m and e > are the electron mass and charge, respectively, iTzli is Planck's constant, jXe = ^jUb/2 is the electron magnetic 
?H moment in terms of jXa, the Bohr magneton, and the electron spin g-factor. Equation ([T]i is coupled to the dynamical Maxwell 
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through the current J = + V x M, which has "free" and spin magnetization parts: 

J/ := fd\d^s\f 
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Extension to muhiple species is straightforward. 

Note, Eqs. ([1]), (|2]), and (O, with (|4| and (|5]l, are to be viewed classically and consequently a full nine-dimensional phase space 
integration, d'^z — d^xd^vd^s, is considered for /. Later we will see how a spin quantization constraint can be applied. 

The Hamiltonian functional for the theory is 



H[^^BJ]=ld\['^.^ + ^-^..B]f+i-Jd\{E'+B^) 
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which can be shown directly to be conserved, but this will become obvious after the Hamiltonian structure is given. 
The noncanonical spin Maxwell- Vlasov bracket is composed of several parts: 

{F,G}mv = jd\f([Ff.,Gf]c 
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with standard partial derivatives denoted by d,. := d /d\ and functional derivatives by Ff := dF/df, etc. Term (|9]l of { , }smv is 
new and accommodates the spin; it is not surprising that it has an inner bracket based on the so(3) algebra ([47]). The remaining 
terms dTji, dHJ, ( fTOl i. and (fTTl i produce the usual Vlasov-Maxwell theory Jim 113, Hzl UMl • It is a simple exercise to show that 
Eqs. ([1}, (l2]l, and (O are given as follows: 
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= {E,H}sMV ■ 



This is facilitated by the identity / d'^zf[g,h] = — f d'^zg[f,h], which works for all three brackets of (fT2l i. ( fTST l. and (fT4l) . 

There are two approaches to obtaining a Hamiltonian description. The usual way is by constructing an action principle by 
postulating a Lagrangian density with the desired observables and symmetery group, and then effecting a Legendre transforma- 
tion, when possible, to obtain a Hamiltonian theory. Alternatively one can postulate an energy functional and Poisson bracket 
as we have done here. When exploring new territory with this latter approach, one must prove directly the Jacobi identity 
{{F,G},H} + {{G,H},F} + {{H,F},G} = for all functionals F, G, and H. With the former approach this is guaranteed if 
the action principle and Legendre transform exist and one can perform a chain rule calculation to obtain a bracket in terms of 
the desired observables. This was done for the Maxwell- Vlasov bracket in |24], where it is necessary to assume the existence of 
a vector potential. However, with the bracket approach one need not assum the existence of a vector potential, in which case the 
Maxwell- Vlasov bracket satisfies 



{{FtG}mv,H}mv 

+ cyc = jd^fV-B 
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This result was quoted in fl^ - see jsTll for a recent recounting of the expHcit (and tedious) details of this early calculation. 
Thus, although the Maxwell- Vlasov Hamiltonian theory is gauge-free, it requires V • B = 0. 

One can construct an action principle for the spin Maxwell- Vlasov theory of the form of ll37l - [39ll and then proceed to the 
bracket {F,G}smv (see e.g. |31, 51]), but we find it easier to prove the Jacobi identity directly. Writing {F,G}smv = {F,G}mv + 
{F, G}s and using :=: to denote the cyclic sum we have 

{{F7G}sMV,H}sMV ■=■ {{F,G}mv,H}mv + {{F,G}s,H}mv + 

+{{F,G}mv,H},, + {{F,GUH}s 
:=: {{F,G}„H}mv + {{F,G}mv,H}s, (16) 

where the second equality follows because of (fTsT i (assuming solenoidal B) and the fact that {F, G}s is a Lie-Poisson bracket 
(see e.g. 1231 129I1 V Thus it only remains to show that the cross terms cancel, which is facilitated by a theorem in (2^; viz., 
when functionally differentiating {F, G}mv and {F, G}s, which are needed when constructing the cross terms, one can ignore 
the second functional derivative terms. These cancel by virtue of the symmetry of the second variation and antisymmetry of the 
bracket. Using the symbol = to denote equivalence modulo the second variation terms, we obtain 

^i^^ ^ [Ff,Gfl. + [Ff,Gf]s + ^-^{FE-d,.Gf^GE-d,Ff) (17) 

^ [Ff.Gfl. (18) 
while all other needed functional derivatives vanish. Thus 
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The first lines of ( fT9] l and ( l20l i cancel by virtue of the Jacobi identities for the brackets [, ]c.b.s on functions, while the second 
line of ( fT9] l cancels upon permutation of the second term. Similarly, the second term of ( |20| ) vanishes. 

Having established the Jacobi identity, we search for Casimir invariants, functionals that commute with all other functionals. 
Using the equations obtained from {C,F} — for all F, we obtain 

Cf' - jd\^{f,s^), (21) 

(22) 
(23) 



= j <i^.x;jQj(x)V-B, 



where JCg, and Kb are arbitrary functions of their arguments. The Casimir is a consequence of the fact that the solution 
to ([T]i is a volume preserving rearrangement, i.e. that the solution can be written as the initial condition on the characteristics. It 
is not difficult to see that ([B can be written in conservation form on the full nine-dimensional space. The dependence of the 
Casimir C^^ is the lift of the so (3) spin Casimir to the kinetic theory. Such inner Casimirs are always give rise to Casimirs of the 
field theory. The Casimir is of course Poisson's equation, an initial condition that would remain preserved should we change 
the Hamiltonian functional. It is a local Casimir because of the arbitrary function K"£(x), which is used here to make the point 
that it is conserved point-wise. The local quantity is technically not the same as the others because its vanishing is required 
for the Jacobi identity. However, this is only technical because {C^,F} ~ 0, for all F, whether or not V • B = 0. 

A consequence of the Casimir Qf^ is that is constant on level sets (contours) of /, which can be viewed as a classical 
prequantization property. If we suppose / has the from / = c(s^)/f(x, v,s,r), then it follows that if fc satisfies ([T]i then / does. 
Choosing 

/ = 5(|s|-/i/2)/,(x,v,s,f) (24) 

we enforce the usual quantization condition and our integrals reduce from integrations over d'^z to d^x,d^vdQ., where dQ. 
denotes the spin sphere as in e.g. l'^. Because of the pure antisymmetry of the so(3) structure constants, Liouville's the- 
orem on characteristics follows immediately; however, for general cosymplectic forms, J, i.e. for brackets of the form 
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[ItS] = d f / dw'J'l(w)dg / dw^ , one can insert a factor of Vdet7 restricted to symplectic leaves to define a proper 'volume' 
measure (see e.g. ll29ll V 

Having found the Casimir invariants we can write down a variational principle for equilibria and then proceed to investigate 
stability by the technique introduce in (l^ (see also |10]), which has become known as the energy-Casimir method (see e.g. 
lll5ll29ll32il] ). First we seek extrema of the quantity := H + C^' + + C^, which must give rise to equations for equilibria; 

S 

— = jr + 47reK-£ + '^/(/,i2)=0 (25) 

— = E-4;rVK-£ = (26) 
5 ^ 

— = B-4;rV^flB = 0, (27) 

where '^f jdf , :— mv^ /2 + 2jj.eS -B/hc and we define the 'particle energy' by (f := J^^ + 4;rfCE. Evidently —4-7:ke 

is the electrostatic potential and B must be an extenal field, i.e. J = V x B = (cf. the results for the Maxwell- Vlasov case 
ll34ll35ll ). Assuming '^^f has an inverse, we obtain the following for the equilibrium distribution function: 

/,(^)='^-i(-^,,2). (28) 

If we chose to be proportional to the usual entropy expression /In/, neglect the dependence on s^, and assume E = 0, an 
acceptible choice, then we obtain the Maxwell-Boltzmann-like equilibrium of |@]. Proceeding to the second variation we obtain 



= -jdh^ff{5ff + — jd\{{8Ef + {8B) 

^ fd\^^ + J- fd'x{{dE)^ + idBf), (29) 



2 J dfe/dS Stt 

where the second equality of ( |29l ) follows upon differentiating the condition + '^^'^ = with respect to /. From ( |29] l we 
immediately draw the formal conclusion that equilibria that are monotonically decreasing functions of S are stable, because 
serves as a Lyapunov functional. More rigorous versions of this have been proved for the Vlasov equation in both the 
plasma and astrophysical contexts (see e.g. ll42ll ). 

Only a limited class of equiUbria come from 5^ = 0, and the reason for this is somewhat subtle. A complete set of equi- 
libria can be gotten from a constrained variational principle with 'dynamically accessible variations' |29], but this will not be 
considered further here. 

In summary, the formulation of an extended kinetic theory for electrons, taking into account the intrinsic spin and the relevant 
magnetization effects, was considered. In particular, the semiclassical limit, valid for length scales large compared to the size 
of the electron wave function was given. Based on the extended phase space, the Hamiltonian structure was discussed, and a 
noncanonical Poisson bracket was found that satisfies the Jacobi identity. Furthermore, we obtained the related Casimir invariants 
and showed the stability of all equilibria with monotonically decreasing distributions. Our findings could act as a guiding tool 
for further extended Hamiltonian theories, including quantum effects from Pauli or Dirac theory, for which a gauge-line has to 
be included in the phase of the definition of the corresponding Wigner function ifTsllssll . Moreover, the stability of the equilibria 
can be an important principle in future numerical studies of strongly magnetized systems. 
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